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ABSTRACT 

Context. Several hundred candidate hybrid pulsators of type A-F have been identified from space-based observations. Their large 
number allows both statistical analyses and detailed investigations of individual stars. This offers the opportunity to study the full 
interior of the genuine hybrids, in which both low-radial-order p- and high-order g-modes are self-excited at the same time. However, 
a few other physical processes can also he responsible for the observed hybrid nature, related to binarity or to surface inhomogeneities. 
The finding that most 6 Scuti stars also show long-period light variations represents a real challenge for theory. 

Aims. We aim at determining the pulsation frequencies of KIC 9533489, to search for regular patterns and spacings among them, and 
to investigate the stability of the frequencies and the amplitudes. An additional goal is to study the serendipitously detected transit 
events: is KIC 9533489 the host star? What are the limitations on the physical parameters of the involved bodies? 

Methods. Fourier analysis of all the available Kepler light curves. Investigation of the frequency and period spacings. Determination 
of the stellar physical parameters from spectroscopic observations. Modelling of the transit events. 

Results. The Fourier analysis of the Kepler light curves revealed 55 significant frequencies clustered into two groups, which are 
separated hy a gap between 15 and 27 d“*. The light variations are dominated by the beating of two dominant frequencies located at 
around 4 d“*. The amplitudes of these two frequencies show a monotonic long-term trend. The frequency spacing analysis revealed 
two possibilities: the puisator is either a highly inclined moderate rotator (v xi 70kms“*, i > 70°) or a fast rotator (v ~ 200kms^*) 
with i 20°. The transit analysis disclosed that the transit events which occur with a sj 197 d period may be caused by a 1.6Rjup body 
orbiting a fainter star, which would be spatially coincident with KIC 9533489. 

Key words. Techniques: photometric - Stars: individual: KIC 9533489 - Stars: oscillations - Stars: variables: delta Scuti - Planets 
and satellites: detection 


1. Introduction 

Recently, ultra-precise photometric data from three space mis¬ 
sions revolutionized the field of variable star studies: the Cana¬ 
dian MO ST satellite (IWa ker et al.ll2003h . the French-led CoRoT 
missi on (Auvergne et al. 2009I) and the NASA mission Ke¬ 


pler (iKoch et al 


•gne et al 

J I2010I) . 


The almost continuous, high-precision 


single-band light curves allow the detection of many pulsation 
frequencies with amplitudes down to the micromagnitude level. 
All these missions provide an enormous amount of extremely 
high-quality data for analysis, enabling state-of-the-art astero- 
seismic studies for many years to come. On the other hand, we 
still need complementary ground-based data for an unambigu¬ 
ous interpretation of the light variations observed by space tele¬ 
scopes, e.g., to determine the physical parameters and modal 
content of the pulsators, or to identify the (mostly non-eclipsing) 
binary/multiple systems. 

A large fraction of intermediate-mass (1.2-2.5 Mq) A and 
F spectral type stars shows pulsations in the intersecting region 
of the main sequence and the classical instability strip in the 


Hertzsprung-Russell diagram (HRD). Among others, we find 
here 5 Scuti (5 Set) variables, pulsating in low-radial-order p- 
and mixed p-g-modes with periods in the 0.01 - 0.3 d range, 
and 7 Doradus (7 Dor) pulsators with high-order g-modes in the 
0.3 -3d period range. The 6 Set pulsations are excited by the 
K mechanism ope rating in the H ell partial ionization zone (see 
e.g. the review bv iBregedbOOOl) . while the 7 Dor pulsations are 
driven by the convective bloc king mechanism o perating at the 
base of the convection zone (IGuzik et alJl2000h . The instabil¬ 
ity regions of these two mechanisms overlap in t he HRD, and 
hybri d 6 Set-y Dor pulsators were first predict ed ( Dupret et al.l 
l2004h . subsequently confirmed dHenrv & Feke 11120^ ). We refer 
to these as hybrids hereafter. The simultaneously excited p- and 
g-modes provide information respectively on the envelope and 
the near core region of the star. Therefore, such targets are very 
important targets for asteroseismic investigations. However, it is 
difficult to find these hybrids with ground-based observations be¬ 
cause the 7 Dor pulsations have usually low amplitudes, and their 
typical frequencies, close to 1-2 d ', are strongly affected by the 
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daily aliasing, as well as by atmospheric effects and instrumental 
drifts. 

Thanks to the recent space-based observations, we now have 
hun dreds of stars classified as A-F-type hyb rid candidates (see 
e.g. llJvtterhoeven et al.ll20fll lHaretCTl2012h . This allows both 
statistical analyses and detailed investigations of individual stars. 
The most striking result of the preliminary analyses of the hy¬ 
brid candidates is that they are not confined to the theoretically 
predicted intersecting region of the two pulsati on types. Instead, 
they occur everywhere in the 6 Set domain (iGrigahcene et al.l 
l2010t[lJvtterhoeven et al.ll201 ih . and it seems tha t practically al l 
5 Set stars show low-frequency light variations dBalonal [201411 . 
In the case of the hotter (5 Set pulsators (Teff > 7500 K), the 
presumed existence of hybrid pulsations challenges the theory 
of stellar pulsations, because we do not expect that the convec¬ 
tive blocki ng mechanism can operate in the thin convective up¬ 
per layer (lBalonall2014ll . Furthermore, most of the candidates 
show frequencies between about 5 and 10d ', even though this 
region has formerly been considered as a ‘frequency gap’, at 
le ast for low spherical de gree (/ = 0-3) modes (see Fig. 2 
of IGrigahcene et aT]l2010tl . This observational finding could be 
explained by the combined effects of frequency shifts of high- 
order g-modes caused by the Coriolis force, det ection of high- 
spher ical-degree modes, and rotational splitting (iBouabid et alJ 
120131) . 

What mechanisms are responsible for the simultaneous oc¬ 
currence of low- and high-frequency variations revealed in so 
many A-F-type stars? Different possible scenarios have to be 
considered; 


1 . Genuine 6 Set-yDor hybrids: both 6 Set and y Dor-type pul¬ 
sations are self-excited in these objects. About a dozen such 
confirmed objects have been investigated in detail, but there 
must be many more genuine hybrids among the known can¬ 
didates. 


2. Eclipsing or ellipsoidal binaries: the object might be a 5 Set 
pulsator in an eclipsing or ellipsoidal binary system. The 
long periodicity in the light variations corresponds to the or¬ 
bital period or half of it, depending on the case. The physical 
parameters of the binary components can be derived inde¬ 
pendently from any theoretical (asteroseismic) assumption. 
These parameters can be used as constraints for subsequent 
asteroseismic investigations. 

3. Tidal excitation: a (5Set star in a close, eccentric bi¬ 
nary, where the external tidal forces excite g-mode pulsa- 
tions (IWillems& AertsllIOO^ . HD 209295 is an example 
dHandler et alJl2002l) . 


4. Stellar rotation coupled with surface inhomogeneities: these 
can also cause detectable signals in the low-frequency re¬ 
gion. Low temperature spots may be an obvious explana¬ 
tion, but mechanisms generating such features are not ex¬ 
pected to operate in normal (non-Ap, non-Am) A-type stars, 
since their convective layers are thin and thus likely have 
weak magnetic fields. N evertheless, s tatistical investigation 
of Kepler A-type stars (lBalonair201 ih shows that the low- 
frequency variability could be caused by rotational modula¬ 
tion in connection with starspots or other corotating struc¬ 
tures. Low temperature spots cause characteristic light-curve 
variations. One low-frequency signal and several of its har¬ 
monics are expected in the Fourier spectrum. However, finite 
starspot lifetimes and/or differential rotation combined with 
spot migration result in a number of significant peaks around 
the rotational frequency on longer time scales. 


However, in the case of KIC 5988140, none of the above 
mentioned simple scenar ios could explain the observations in 
a fully satisfactory way (iLampens et al.ll207^ . Indeed, in this 
case, the scenario of binarity has been discarded on the basis of 
the double-wave pattern of the radial velocity curve mimicking 
the pattern of the light curve (after prewhitening of the higher 
frequencies). The other two remaining scenarios explaining the 
low frequencies in terms of g'-mode pulsation or rotational mod¬ 
ulation are also problematic since the predicted light-to-velocity 
amplitude ratio based on two simple spotted surface models is 
40 times larger than observed. 

KIC 9533489 is another Keple r hybrid candidate investi¬ 
gated by lUvtterhoeven et aP (1201 ill , and which turns out as an 
exotic case: in addition to its confirmed hybrid character, we also 
detected a few transit events in the Kepler light curve. In the fol¬ 
lowing sections, we present the results of the frequency analyses 
and derived (physical and pulsational) properties of this object, 
as well as a possible model for the transit events detected in the 
light curve. 


2. The Kepler data 

KIC9533489 (Kp = 12.96 mag, 02000 = 19*'38”4L7^ ^2000 = 
-i-46‘^07™2L6’‘) was observed both in long-cadence (LC) and 
short-cadence (SC ) mode. The LC dat a have a sampling time 
of 29.4 minutes (IJenkins et al.l l2010l) . while the SC obser- 
vation s provide a brightne ss measurement every 58.85 sec¬ 
onds (iGilliland et al.l l2010l) . We used the ‘raw’ flux data of 
KIC 9533489 available from the Kepler Asteroseismic Science 
Operations Center (KASOC) database. We analysed the LC data 
of 17 quarters of Kepler observations (Q1-Q17, 13th May 2009 
- 11th May 2013) and the SC data of 3 quarters (Ql, Q5.1-3, 
Q6.1-3, 13th May 2009 - 22nd September 2010). The contami¬ 
nation level of the measurements is estimated to be only 1.6 per 
cent. 

The strong long-term instrumental trends present in the light 
curves were removed before analysis. First, we divided the time 
strings into segments. These segments contain gaps no longer 
than 0.1-0.5 d. Then, we separately fitted and subtracted cubic 
splines from each segment. We used one knot point for every 
1000 (SC data) or 100-200 points (LC data) to define the splines. 
We omitted the obvious outliers and the data points affected by 
the transit events. The resulting SC and LC light curves consist 
of 307950 and 62418 data points, spanning 498 (with a large 
gap) and 1450 days, respectively. 


3. Spectroscopy 

For an unambiguous interpretation of the results of the light- 
curve analysis, the atmospheric stellar parameters, i.e. the ef¬ 
fective temperature (Tefr), surface gravity (logg) and projected 
rotational velocity {v sin /), must be known. Therefore, we ob¬ 
served KIC 9533489 in the 2013 and 2014 observing seasons 
with the High Effi ciency and Resolut ion Mercator Echelle Spec¬ 
trograph ( HERMES. [Raskin et al.l201 ill attached to the L2-m Mer¬ 
cator Telescope, with the Elbre-fed Echelle Spectrograph (eies) 
at the 2.5-m Nordic Opt ical Telescope (NO T), and also with 
the INTEGRAL instrument (lArribas et al.l 1 199^ using the 4.2-m 
William Herschel Telescope (WHT). All instruments are lo¬ 
cated at the Roque de los Muchachos Observatory (ORM, La 
Palma, Spain). We also acq uired data with th e High Resolution 
Echelle Spectrometer (hires. [^gt et al.ll994t) attached to the 10- 
m KeckI telescope (W. M. Keck Observatory, Hawaii, U.S.). Ta- 
ble[T]shows the journal of the spectroscopic observations. 
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Table 1. Log of spectroscopic observations of KIC 9533489. 


Telescope 

Instrument 

Wavelength range 
(nm) 

Resolution 

Obs. Time 
(JD) 

Exp. time 

(s) 

Obs. 

# 

S/N 

1.2-m Mercator 

HERMES 

380-900 

85 000 

2456514.4 

1800 

3 

20 

2.5-m NOT 

FIES 

370-730 

45 000 

2456574.4 

1800 

3 

55 

4.2-m WHT 

mXEGRAL 

570-700 

2750 

2456783.7 

600 

1 

70 

10-m KeckI 

HIRES 

386-678 

67000 

2456 845.0 

1132 

1 

100 


Notes. S/N is the signal-to-noise ratio per sample of the mean spectrum computed in the A = 496 - 564 nm wave length range (hermes, fies and 
HIRES spectra), or at /t ~ 570 nm (integral spectrum), and by adopting a procedure described in IStoehr et alj l l2007h . The hires spectrum has gaps 
in the 477^98 and 628-655 nm domains. 


We obtained the reduced spectra from the dedicated pipelines 
of the HERMES, FIES and hires instruments, while the single-order 
WHT spectrum centered on Ha was reduced using iraiQ rou¬ 
tines. All spectra were corrected for barycentric motion, then 
combined into a single averaged spectrum per instrument. 

Though the signal-to-noise ratio (S/N) of the averaged her¬ 
mes spectrum is low (S/N = 20 per wavelength bin at A ~ 
500 nm), calculating a cross-correlation function with an FO- 
star synthetic spectrum, excluding regions of B aimer and tel¬ 
luric lines, we were able to derive the radial velocity (RV) for 
that epoch. We performed the same procedure on the fies and 
HIRES spectra, too. We found that the RV measurements are con¬ 
sistent with each other within the errors, pointing to a value of 
6.5 + 1.5kms That is, the star showed no sign of RV change 
over a period of ~ 330 days, and therefore cannot be considered 
as a multiple system. 

Because it has a sufficiently high S/N and that it covers a 
wide wavelength domain that contains all the majour Balmer 
lines, we started to derive the star’s astrophysical parameters 
on the FIES averaged spectrum. We estimated vsini by adopt- 
ing the Fourier transfor m (FT) approach (for an overview, see 
iGravll 1 99^ lRoveiil2005L and references therein). To increase the 
S/N and to limit the impact of blends, the FT was performed on 
an average line profile obtained through cross-correlation of the 
spectrum with an FO mask from which we excluded the strongest 
lines. The result shows that the star is probably a moderate rota¬ 
tor (see TableO. 

We estimated by fitting the Ha, H/1, Hy, and Hb line 
profiles of the fies spectrum. The continuum normalisation of 
the spectrum was realised in three steps. A first normalisation 
was done with the iraf/continuum task to perform a very prelim¬ 
inary determination of Teff and logg in the 400^50 nm range 
(metallicity was kept Solar and microturbulence was fixed to 
2 km s“'). The fit was realised by using the girfit program based 
on the minuit minimization package and by interpolating in a 
grid of synthetic spectra. These synthetic spectra and corre- 
spon ding LTE model atmos pheres were computed with the syn- 
SPEC dHubenv & Lanzil 19951 an d references therein), and with at- 
las9 (ICastelli & Kuruc3 20041 and references therein) computer 
codes, respectively. They were then convolved with the rota¬ 
tional profile and with a Gaussian instrument profile to account 
for the resolution of the spectrograph. The observations were 
then divided by the closest, not normalised, synthetic spectrum. 

After being sigma clipped to remove most mismatch fea¬ 
tures, the resulting curve was divided into 2 to 4 smaller re¬ 
gions which were smoothened with polynomials of degrees 5 


* iRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 


Table 2. Effective temperature values derived from the fit of the hydro¬ 
gen line profiles observed with different instruments. 


Instrument 

Teff (K) 

Line profiles considered 

HERMES 

7250 + 150 

Ha, H/3, Hy 

FIES 

7350 + 200 

Ha, H/3, Hy, H6 

INTEGRAL 

7200 + 20 

Ha 

HIRES 

7190 + 70 

HS,Hy 


to 10. The observe d spectrum was divi ded by the resulting func¬ 
tion (see Fig. 1 of iFremat et al.l (l2007h for an illustration of the 
procedure) then normalised with the iraf/continuum task and a 
low order spline. Regions around Ha, li/3, Hy, and H5 were then 
isolated and fitted separately. Because below 8000 K, hydrogen 
lines are only sensitive to effective temperature, all other param¬ 
eters were kept fixed and we tried as much as possible to exclude 
from the fit most metal lines. Each time, 3 to 5 different starting 
points were chosen between 6000 and 8000 K, which resulted in 
a median value and an inter-quantile dispersion of 7350 + 200 K. 

The same procedure was then applied on the other spectra, 
except that fewer hydrogen lines were available. We list the val¬ 
ues we obtained in Table |2l where the last column details the 
hydrogen lines that we considered. For the hermes and fies spec¬ 
tra the error bars were deduced from the scatter of our different 
determinations and trials, while for the integral and hires deter¬ 
minations we adopted the numerical value provided by the mi¬ 
nuit package when it computes the covariance matrix. The final 
Teff value is obtained by combining these different estimated into 
a rounded mean which error is provided by the largest deviation 
(see TableO. 

Since our initial scope was not to perform a detailed chemi¬ 
cal abundance analysis, we derived the metal content by scaling 
the abundance of all elements heavier than boron by the same 
factor which we call here after ‘metallicity’. For this purpose we 
studied the region between 496 and 564 nm. To identify the fea¬ 
tures that are the most sensitive to surface gravity, we considered 
two synthetic spectra of same effective temperature (7200 K) and 
V sin i but having different extreme log g values (3.5 and 5.0). The 
difference between the two spectra exceeds the limits given by 
the S/N of the observed spectra only in the region around the 
Mg I triplet (515 to 520 nm). We then constructed a first regu¬ 
lar grid of synthetic spectra for different values of the microtur- 
bulent velocity and metallicity. The xi value between the ob¬ 
servations and the spectra of each (timic, metallicity) node was 
computed from 496 to 564 nm, excluding the 515-520 nm part. 
Finally, a Levenberg-Marquardt minimization algorithm was ap¬ 
plied around the minimum and a covariance matrix was com¬ 
puted to estimate the error bars by interpolating in a second grid 
of spectra having smaller steps in metallicity and microturbu¬ 
lence. Log g was then derived by fixing the microturbulence and 
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Fig. 1. Observed hires spectrum (in black) and synthetic model (in red) 
of an FO-type star with usini = 70km s“' over the wavelength range 
400-455 nm. The residual spectrum is drawn in blue. 


metallicity to the best values found in the previous step, and by 
applying the Levenberg-Marquardt algorithm to the region be¬ 
tween 515-520 nm. Because, the final log^ value did not have 
any significant impact on the agreement between synthetic spec¬ 
trum and observations in the other considered wavelength do¬ 
mains we stopped the process. 

Part of the hires spectrum is shown in Fig. [1] where it is com¬ 
pared to the synthetic spectrum that results from our atmospheric 
parameter determination procedure. The residuals between the 
two spectra are also plotted in the same figure. These residu¬ 
als are significant compared to the S/N. We checked therefore 
whether we could find any trace of an additional contribution in 
the spectrum as suggested by the transit modelling (see Sect. 

First, we computed the cross-correlation of the fies spectrum 
with a G5 spectral template. We did not find any evidence of 
absorption lines due to a cooler companion. This could mean that 
the presumed companion’s contribution to the total light is too 
low to be detected in the spectrum and/or that its spectral lines 
are too broadened (i.e. due to rotation or/and macroturbulence). 
A light factor corresponding to Am = 2.25 mag might explain 
such non-detection. 

We next computed an (F0-I-G5) composite model spectrum 
and looked at the variance between this model and the hires 
spectrum in the wavelength range 500-560 nm for a wide range 
of possible radial velocities. The G5-synthetic spectrum was 
computed assuming ttsint = lOkms The minimum variance 
was found around 31kms"', i.e. different from the radial ve¬ 
locity of the Kepler star. This second test might hint toward the 
existence of an additional stellar object, but we remain unable to 
identify the lines of the transit host directly in the spectrum. 

Table [3 summarizes the stellar parameters derived from the 
spectra. Knowing the values of and log g, and assuming a so¬ 
lar metallicity, we estimated the mass, luminosity, radius and age 
of KIC 9533489 (see also in Table [T|l app l ying th e evolutionary 
tracks and isochrones from iGirardi et al.l (|2000|) . The Teir and 
logg values indicate that KIC 9533489 is situated in the over¬ 
lapping region of the 6 Set and y Dor observational instability 


Table 3. Stellar parameters of KIC 9533489, determined from spectro¬ 
scopic ob servations (‘spectra ’) and evolutionary tracks and isochrones 
(‘models’, IGirardi et alj200d) . 



SPECTRA 


MODELS 

Teff 

7250+ 100 K 



log0 

4.25 ± 0.2 dex 


5.8!f^Lo 

V sin i 

70±7kms-' 

Rt 

1.53:“^ Ro 

^mic 

2.9 + 0.3kms“' 

Age 

350!™ Myr 

[M/H] 

-0.05 + 0.2 

P* 

0.42!0;i>o 

Kad 

6.5 + 1.5kms“^ 





Fig. 2. Location of KIC 9533489 on the surface gravity vs. effective 
temperature d iagram. Evolutionar y tracks for 1.4, 1.5, 1.6, and 1.7 Mq, 
as derived bv IGirardi et all ( l200Qll . are plotted with solid lines. Dotted 
lines show 63 Myr-1.12 Gyr isochrones. 


strips (see e.g. fig. 10 of lUvtterhoeven et al.ll201 ll) . The age esti¬ 
mation suggests a young object, located near the zero-age main 
sequence (ZAMS), but the uncertainty is very large (Fig.|2l. 


4. Adaptive optics observations 

We obtained near-infrared adaptive optics (AO) images of 
KIC 9533489 in the Br-y (2.157 pm) filter using the NIRC2 im¬ 
ager behind the natural guide star adaptive optics system at the 
Keckll telescope on 04 September 2014. Data were acquired 
in a three-point dither pattern to avoid the lower left corner of 
the NIRC2 array which is noisier than the other three quadrants. 
The dither pattern was performed three times with step sizes of 
3", 2.5", and 2.0" for a total of nine frames; individual frames 
had an integration time of 30 s. The individual frames were flat- 
fielded, sky-subtracted (sky frames were constructed from a me¬ 
dian average of the 9 individual source frames), and co-added. 

The average AO-corrected seeing was FWFIM = 0.048"; the 
pixel scale of the NIRC2 camera is lOmas/pixel. 5cr sensitiv¬ 
ity limits were determined for concentric annuli surrounding the 
primary targets; the annuli were stepped in integer multiples of 
the FWHM with widths of 1 FWHM. The sensitivity achieved 
reached Amag » 8 mag at 0.2" separation from the central tar¬ 
get star in the ba nd of observation. U sing typical stellar Kepler- 
infrared colours (iHowell et al.ll201^ . this value translates into 
the Kepler bandpass to be about Amag a; 10-11 mag. 

One source was detected approximately 1.1" to the west of 
the primary target (Acr = 1.12 + 0.005"; A6 - 0.09 + 0.005"; 
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Fig. 3. Keck II adaptive optics image of KIC 9533489. 
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Fig. 0. The companion is 3.412 + 0.014mag fainter than the 
target at 2.2jum; deblending the 2MASS Ks photometry, we 
hnd the infrared brightnesses of the target and the companion 
to be Ks - 12.055 + 0.020mag and Ks = 15.467 + 0.024mag. 
No observations in other wavelengths wer e obtained, but base d 
upon the typical Kep/er-infrared colours (iHowell et al.ll2012|) , 
the companion has a Kepler bandpass magnitude of Kp - 17.9 + 
0.8 mag; thus, the companion is > 5 magnitudes fainter than the 
primary target and contributes less than 1% to the Kepler light 
curve. 


5. Frequency analysis 

We performed the Fourier analysis of both the SC and LC data 
sets and compared the frequencies obtained. The increased sam¬ 
pling of the SC data allows the direct determination of pulsa¬ 
tion frequencies above the LC Nyquist limit (/uc.Nyq = 24.5 d"\ 
/sc.Nyq = 734d“'). The LC observations, having a longer time 
base, provide better frequency resolution and allow the inves¬ 
tigation of the long-term frequency stability. The Rayleigh fre¬ 
quency resolutions of the long- and short-cadence data sets are 
/lc.r = 0.0007 and /sc.r = 0.002 d“', respectively. 

We performed a standard successive prewhitening of 
the light curves fo r signihcant signals usin g Period04 
dLenz & Breged[^05h and SigSpec (lReegenll2007h . The results 
were also checked with the photometry modules of famias dZimal 
I2OO8I) . The frequencies obtained were consistent within the er¬ 
rors using the different programs. Fig. |4] presents the Fourier 
transforms of the original and the pre-whitened SC data sets. We 
stopped the process whe n none of the rema ining peaks reached 
the threshold of S/N « 4 dBreger et alJll993h . Altogether, 63 fre¬ 
quencies were determined in this way. We also checked the spec¬ 
tral significances (sig) of the frequencies, as derived by SigSpec. 
All of the peaks were found to be above t he default thre shold 
value for determining frequencies (5.46, see lReegenll2007h . The 
lowest amplitude peaks have sig~8 (cf. Table |4]i. We checked 
which of these frequencies can also be found in the LC data. Be¬ 
low the LC Nyquist limit, we found all but eight frequencies. All 
of the eight missing frequencies are below 1 d ~^ Both the in- 
strumental effects in the original data (see e.g. lMurDhvll2Q12h 


Fig. 4. Fourier transforms of the original SC light curve and of the light 
curve pre-whitened for 2, 9, 46, and 63 frequencies. The blue dashed 
lines denote the 4(A) significance level, where (A) is calculated as the 
moving average of radius 8d“' of the pre-whitened spectrum. The fre¬ 
quency gap between x:15 and 27 d“* is indicated with gray band. The 
vertical dotted line denotes the Nyquist limit of the LC data. 


and the correction method we applied to remove these affect 
this low-frequency region. Therefore, we excluded these eight 
peaks from the list of accepted frequencies. Checking the re¬ 
maining significant frequencies using the Kepler Data Charac¬ 
teristics Handbook (KSCI-19040-004), none of them was found 
to be of instrumental origin. 

Note that the frequencies detected in the SC data above the 
LC Nyquist limit, except for the three lowest amplitude ones, 
also appear in the LC data. These can be identihed as Nyquist 
aliases at fi t ,r = 2 fi ,r Nvn - /i,sc. split by Kepler's orbital fre¬ 
quency (iMumhv et al.ll2013l) . Fig. |5]presents the FT of the origi¬ 
nal SC and LC data, respectively. We can detect both real pulsa¬ 
tion frequencies and such Nyquist aliases above lL3d“' in the 
Fourier transform of the LC light curve. 

Table |4] lists the 55 accepted frequencies and the differences 
between the frequency values derived from the SC and LC data, 
too. They differ only in the fourth-sixth decimals. Note that the 
frequency difference (0.0056 d“^) of the closest peaks (/m and 
/51) is 2.8 times larger than the Rayleigh frequency resolution 
of the SC data set, so all of the identihed frequencies are well- 
resolved. 

Two closely spaced peaks at 4d“' (/i and /2) dominate 
the light variation of KIC 9533489 (see Fig.0. Most of the peaks 
(42) are below 15 d“ *, and 27 of these are below 5 d“'. The rest of 
the frequencies (13) are in the 27-38 d ' range. There is a clear 
frequency gap between 15 and 27d '. We indicate this region 
with a gray band in Figs. |4] and |5] Such gaps play an important 
role in hybrid star studies, as these allow to distinguish between 
different mode-driving mechanisms. 

We checked the frequencies for possible combination terms 
(k/i or pf\ + qf\). Co nsidering test results on simulated light 
curves (iPanicsl 1201^ . we restricted our search only to the 
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Fig. 5. Fourier transforms of the original SC (upper panel) and LC (bottom panel) light curve, respectively. The frequency gap between x:15 and 
27 d“* found in the SC data is indicated with gray band. The vertical dotted line denotes the Nyquist limit of the LC data. Red lines connect the 
lowest and the highest frequency peaks of the SC data with their Nyquist aliases detected in the LC Fourier transform. 
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Fig. 6. A section of the SC light curve of KIC 9533489 and the fit of the 
two dominant frequencies (red line). The plot exhibits the beating effect 
dominating the light variations. 


second-order combinations, i.e. \p\ - 1 and |q'| = 1. We accepted 
a peak as a combination if the amplitudes of both parent frequen¬ 
cies were larger than that of their presumed combination term, 
and the difference between the observed and the predicted fre¬ 
quency was not larger than the Rayleigh resolution of the SC data 
(0.002 d“'). This way, we found one harmonic peak (/49 = 2/2) 
and 11 other combination terms. Note that although /30 is close 
to 2/1, it does not fulfil the above criterion. We indicate the com¬ 
bination terms in Table H] Beside the combinations consisting 
of independent parent frequencies, we list two additional cases 
in parentheses in Table 01 too. These are equal to two and four 
times the value of f 2 s, respectively, which is itself a combination 
of two large-amplitude frequencies (/28 - fi— fi)- It is not clear 
yet, whether this is a coincidence, or whether f 2 s (0.84d“') has 
a particular importance for some unknown reason. 

5.1. Stability investigations 

Strong amplitude variability on month- and year-long timescales 
was reported in the case of the A-F-type hybrid KIC 8054146 
(iBreger et al.ll2012l) . We also examined the long-term stability 
of the most dominant peaks. We used the LC data set, because 
it spans a much longer time interval than the SC data. We in¬ 
vestigated the nine highest-amplitude frequencies below the LC 
Nyquist limit. First, we calculated the actual frequencies and 
their amplitudes for each quarter (Q1-Q17). Then, we deter¬ 
mined values that characterise the observed variability of 
these parameters as the deviation of the quarterly data points 
from their time-averages, relative to the uncertainties of the mea¬ 
surements. We did not detect variability in any of the nine fre¬ 
quencies, as <1.3 for each of them. 

However, this test revealed the variability of two amplitudes. 
We obtained values for the amplitudes of the two most dom- 
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Fig. 7. Amplitudes of /i and /2 in the different light-curve segments of 
the Q1-Q17 LC data. The dashed line shows the error-weighted average 
of the data points, the continuous line is a linear fit to the data. 

inant frequencies, f\ and /2, of 8.3 and 4.8, respectively. The;^f^ , 
value of the other seven amplitudes remained below 1.5. Fig.lT] 
presents the amplitude values of fi and /2 for Q1-Q17. The am¬ 
plitudes show a linear trend and vary similarly but in opposite 
directions. The continuous lines are linear fits to the data points, 
which show a 2.4 and 2.0 per cent relative amplitude variation 
per year. values for the linear fits are 1.1 and 0.6, for 

/i and /2, respectively. These amplitude changes could be ex¬ 
plained in terms of energy exchange between the two dominant 
modes and interpreted as coupling. 

We also checked whether we can detect any sign of phase 
modulation ( PM) of the dominant modes, applying the method 
presented by iMumhv et al.l (l2014l) . A periodic phase modula¬ 
tion would indicate the presence of a stellar or even a substel- 
lar companion of KIC 9533489. The panels of Fig. [8] show the 
time delays of /i and /2, calculated from 10-days segments of 
the LC data, and their Fourier transforms, respectively. The for¬ 
mula used to calculate the time delays is Tij = ((pij- (f)j)l{27TVj), 
where (pj is the mean of the (pj{t) phases derived f rom the differ¬ 
ent lig ht curve segments for a fixed Vj frequency dMurnhy et al.l 
I2OI4I) . None of these point to any periodic variation, that is, 
KIC 9533489 is not a member of a 20 d < Porb < 4 yr binary sys¬ 
tem: neither the time delay plots, nor the corresponding Fourier 
transforms have a common significant peak corresponding to a 
possible common orbital period. 

5.2. Frequency and period spacings 

In the case of a genuine hybrid, KIC 9533489 will show simulta¬ 
neously excited high-order g- and low-order p-modes. Thus, we 
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Table 4. Frequency content of the KIC 9533489 Kepler data. 



Frequency 

[d-i] Ar 

[pm 

Ampl. 

[PPt] 

Phase 

[27t] 

S/N 

sig 

Combinations 

/l 

4.007089(4) 

-0.04 

46.37834 

0.333 

0.566 

186.2 

9967.3 


fl 

4.443236(4) 

0.10 

51.42635 

0.322 

0.039 

180.3 

10647.4 


h 

3.60269(1) 

-0.23 

41.6978 

0.120 

0.185 

67.3 

1747.5 


/4 

2.87088(1) 

-0.33 

33.2277 

0.116 

0.947 

64.6 

1693.5 


h 

3.68992(1) 

0.12 

42.7075 

0.113 

0.940 

63.1 

1646.5 


fe 

10.08357(1) 

0.27 

116.7080 

0.102 

0.850 

59.0 

1385.9 


fl 

32.58125(2) 


377.0979 

0.080 

0.506 

46.7 

850.6 


h 

31.65313(2) 


366.3557 

0.065 

0.275 

38.0 

580.5 


/9 

30.19377(2) 


349.4650 

0.064 

0.562 

37.2 

563.3 


/lO 

3.77254(2) 

0.04 

43.6637 

0.057 

0.884 

31.6 

446.9 


/ll 

32.35529(2) 


374.4825 

0.055 

0.554 

32.4 

432.9 


/l2 

6.36563(2) 

0.18 

73.6763 

0.053 

0.597 

29.8 

395.0 


/l3 

1.15926(3) 

-0.05 

13.4173 

0.052 

0.243 

28.7 

374.6 


/l4 

1.87679(3) 

-0.22 

21.7221 

0.044 

0.401 

24.7 

280.3 


/l5 

2.50313(3) 

0.38 

28.9714 

0.042 

0.670 

23.4 

253.5 


/l6 

1.51617(5) 

1.11 

17.5483 

0.041 

0.103 

23.0 

254.1 


/l7 

29.73612(3) 


344.1680 

0.039 

0.163 

22.7 

217.2 


/l8 

8.75884(4) 

0.71 

101.3754 

0.034 

0.383 

19.6 

173.6 


/l9 

3.38343(5) 

-0.99 

39.1601 

0.029 

0.698 

16.0 

121.7 


fio 

7.37521(5) 

-0.66 

85.3612 

0.028 

0.557 

15.8 

115.7 

fl + /lO 

fll 

35.37792(5) 


409.4667 

0.025 

0.664 

14.8 

94.6 


fll 

37.60866(5) 


435.2854 

0.024 

0.258 

14.3 

86.0 


/23 

11.92258(6) 

0.15 

137.9929 

0.022 

0.674 

13.0 

72.6 


/24 

8.57253(6) 

0.95 

99.2191 

0.022 

0.654 

12.3 

70.2 


/25 

27.04778(6) 


313.0530 

0.021 

0.033 

12.1 

62.9 


fie 

33.05026(7) 


382.5261 

0.019 

0.553 

11.0 

51.5 


fn 

3.89026(7) 

-0.58 

45.0261 

0.019 

0.298 

10.5 

51.2 


fis 

0.83999(7) 

2.10 

9.7221 

0.019 

0.780 

10.2 

49.5 

fi-fi 

f29 

1.67951(7) 

-1.71 

19.4388 

0.018 

0.575 

9.8 

46.1 

(2/28) 

fio 

8.02468(7) 

0.40 

92.8782 

0.017 

0.135 

9.7 

43.6 


/ 3 I 

12.39149(8) 

-1.47 

143.4200 

0.016 

0.053 

9.5 

39.1 


fn 

5.22971(9) 

-0.04 

60.5291 

0.015 

0.166 

8.4 

33.3 


fn 

5.60943(9) 

1.87 

64.9240 

0.014 

0.787 

7.9 

29.4 


fn 

1.56217(9) 

-1.43 

18.0806 

0.014 

0.512 

7.9 

28.4 


fn 

14.5764(1) 

-0.37 

168.708 

0.012 

0.650 

7.2 

22.6 


fie 

3.6622(1) 

-0.95 

42.387 

0.012 

0.135 

6.9 

22.1 

hi + /i5 

fn 

33.6440(1) 


389.398 

0.012 

0.977 

6.9 

19.9 


fn 

2.5311(1) 

-0.72 

29.295 

0.010 

0.532 

5.7 

15.9 

h - hi 

fn 

7.9278(1) 

-1.43 

91.757 

0.010 

0.868 

5.8 

16.2 

fn + fn 

/40 

6.6502(1) 

2.25 

76.970 

0.010 

0.425 

5.8 

15.5 


/41 

7.7748(1) 

-1.34 

89.986 

0.010 

0.261 

5.6 

14.7 


/42 

35.0501(1) 


405.673 

0.010 

0.286 

5.8 

14.4 


fn 

1.1769(2) 

0.69 

13.622 

0.009 

0.357 

5.1 

13.1 


/44 

2.1616(1) 

0.65 

25.018 

0.009 

0.397 

5.1 

12.7 

fn - h 

/45 

4.8035(1) 

-0.17 

55.596 

0.009 

0.145 

5.0 

11.9 

hi - fn 

/46 

1.2225(2) 

-1.07 

14.150 

0.009 

0.514 

4.8 

10.9 

hi - h 

/47 

3.3593(2) 

0.75 

38.881 

0.008 

0.586 

4.6 

9.9 

(4/28) 

fn 

1.1976(2) 

4.00 

13.861 

0.008 

0.381 

4.5 

8.0 


fn 

8.8862(2) 

-2.26 

102.850 

0.008 

0.775 

4.7 

9.8 

2fi 

/50 

34.1421(2) 


395.163 

0.008 

0.003 

4.8 

9.9 

h + fn 

/51 

1.5106(2) 

-5.66 

17.483 

0.008 

0.180 

4.5 

8.6 

fe - fi4 

hi 

2.5667(2) 

0.99 

29.707 

0.008 

0.992 

4.4 

9.2 

fl - /l4 

hi 

2.2978(2) 

3.10 

26.595 

0.007 

0.382 

4.1 

8.0 


/54 

31.8377(2) 


368.492 

0.007 

0.223 

4.3 

8.0 


fn 

1.2767(2) 

0.03 

14.776 

0.007 

0.494 

4.1 

7.8 



Notes. The errors given in parentheses are standard uncertainties derived from the least-squares fitting. We present the frequencies, amplitudes 
and phases derived from the SC light curve, as the SC data allow the direct detection of peaks above 24.5 d“*. Furthermore, we list the freq uency 
differences of the SC and LC frequencies in the third column. The spectral significances (sig) were calculated using SigSpec (lReegenll2007ll . 

<‘'>Ay; = (y;,sc-/:,Lc)xiond-‘]. 
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Fig. 8. Time delays calculated for the two dominant modes fi and /2 
(top panels), and their Fourier transform (bottom panels). 


checked, whether there are any regularities between the periods 
and frequencies in the presumed g- and p-mode regions, respec¬ 
tively, predicted by stellar pulsation theory. For referenc es on the 
theore ti cal background of non- radial pulsa t ions, s ee e.g. lTasst^ 
(ll980ll . lTrnno^t^(ll989ll and lAerts 61311 ( 1201 ( 111 . 

High-radial-order gr-modes are excited in 7 Dor pulsators. 
According to the first-order asymptotic approximation of non- 
radial stellar pulsations, the separation of consecutive overtones 
having the same { spherical harmonic index (AP( - Pn+i,e-Pn.e) 
is: 


APe - Ho/ y/e{{+l). 


( 1 ) 


where 


Ho = 27t^ 



1 


and N is the Brunt-Vaisala frequency. That is, period spac- 
ings are expected to be constant in a chemically homogeneous 
star, but inhomogeneities (rapid changes of the mean molecular 
weight in the transi tion zones) cause devi ations from this con - 
stant value (see e.g. iMiglio et al.ll2008l an d iBouabid et al.ll 20 Tll) . 
Following Eq. [T] the period spacing ratios for the different { 
modes are: 


AP(=i/APr =2 ^ V3 = 1.73, (2a) 

- V2= 1.41. (2b) 

In the case of low-radial-order p-modes, even though the 
pulsation modes are in the non-asymptotic regime for most of 
the dSct stars, we still can notice a clustering of frequencies 
with the presumed large separ ati on, Av^ = Vn+i.r - V ii.r (see e.g. 
Garcia Hernandez et al.il200^ and lPanaro et alj|2013h . The large 
separation can be calculated in a first approximation as: 

Av 
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where Cj is the sound speed. The equation given for the large sep¬ 
aration (integral over the inverse sound speed) is valid for high 
radial order p-modes (lTassoulll980HGoughll990h . It was gener¬ 
alized in terms of the mean density of the stars bv ISuarez et al.l 
(l2014h . based on numerical calculations using a huge model grid. 
The relation obtained for the 6 Set stars domain can be approxi¬ 
mated by 

Av/Ave = (3) 

where Avq = 134.8 pHz and ^ = 1.408 gcm“^. 

Another well-known relation in the case of 6 Set pulsations is 
the frequency-ra tio of the radial fun damental and first-overtone 
radial mode (see lStellingwerflfl979h : 


/fu//io^ 0.77-0.78. (4) 

In the case of high-overtone (n » €) modes and slow rota¬ 
tion, the frequency differences of the rotationally split compo¬ 
nents can be calculated by the asymptotic relation: 


SVn 4 ,m = Sm{\ - (5) 

where the coefficient C,,/ l/^(^ + 1) for g'-modes, C„/ 0 for 

p-modes, and Q is the (uniform) rotation frequency. As 5v ~ 
[1 - ll({( - 1 - 1 )], the ratio of the rotationally split { - I and { - 2 
g-modes is in the case of equal 6m: 


5vt=\,gl6vt=l,g ^ 0.60. ( 6 ) 

We also have in a first approximation: 


6vp ^ 26vc=\^g. 


(7) 


We investigated the frequency set of KIC 9533489 considering 
these relations. 

Note that in the case of rapidly rotating stars, additional 
terms to Eq. |5] have to be considered in the calculation of ro¬ 
tationally split frequencies, considering higher-order rotational 
effects. In this case, 6vn,t,m also depend s on v„j and on m 
(m^-rn ^). E or a theoretical backgr ound, see lGough & Thompson! 


( I^90|) and Kjeld sen et alJ ( 1998h f or the acous t ic mod es, and 
IChlebow^ (1 19781) and iDziembowski & Goode! (Il992h for the 
g-modes. In Sect. 15.2.21 we use the simple form of Eq. |5]to de¬ 
rive a first approximation for the rotation period of the star. 


5.2.1. Frequency ratios 

Considering the modelling results of ISuarez et al.l (l2014t) on 
main-sequence 6 Set stars, knowing the mean density of 
KIC 9533489, we estimated the theoretical radial fundamental 
mode’s frequency. Taking into account the estimated radius and 
mass, and their uncertainties (see Table O, the mean density of 
the star normalized to the solar value may be between 0.2 and 
0.6 p/pq. These values indicate that the radial fundamental mode 
could lie between ~ 13 - 23 d * (150-260pHz), practically in 
the ob served frequency gap, according to fig. 2 of ISuarez et al.l 
(l2014h . This means that the radial fundamental mode is either 
/as (14.58 d“*), the only mode between 13 - 23 d“', or the radial 
modes are not exceeding the S/N= 4 significance level. There is 
no first overtone frequency pair for f^s with frequency ratio in 
the range of 0.77-0.78 according to Eq.|4] 
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As the radial fundamental mode represents the low frequency 
limit for p-modes, we can assume that the frequencies below the 
gap (up to 14.6d ') are ‘yDor-type’ p-modes. Rapid rotation 
may explain both the shifts of gr-modes int o higher frequencies 
than the ‘classical’ y Dor pulsation domain (iBouabid et al.IZOOh 
and the lack of radial modes. Indeed, rapidly rotating 6 Set stars 
generally tend to show mainly nonradial modes o f lower ampli - 
tude in comparison with the slow rotators (see e.g. lBregerf2000ll . 
It is also possible that some of the observed frequencies above 
5 d ' are non-zonal (m + 0) modes generated by rotational split¬ 
tings. 


5.2.2. Frequency spacings and stellar rotation 

We studied the distribution of the frequencies applying differ- 
ent methods: by ca lculating their Fourier transforms (see e.g. 
[Handler et alJ Il997h. generating his tograms of the frequency 
differences (e.g. in [Matthe^ IZOOTj), and using Kolmogorov- 
Smirnov (K-S) tests (iKawalei Il988h We defined and analysed 
several sets of frequencies: those below 5, 10, 15, and above 
27 d^'. In the Fourier tests, we assigned equal amplitude to all 
frequencies. Local maxima in the FTs indicate that characteris¬ 
tic spacings exist between the frequencies. In the K-S test, the 
quantity Q is defined as the probability that the observed fre¬ 
quencies are randomly distributed. Thus, any characteristic fre¬ 
quency spacing in the frequency spectrum should appear as a lo¬ 
cal minimum in Q. We present the results of these computations 
in Fig. 121 

Provided that both g- and p-modes occur, we expect to find 
periodicities that belong to rotationally split frequencies, accord¬ 
ing to Eqs. |5j | 6 ] and |7] However, note that these relations are 
valid for slow rotators only, and in the case of KIC 9533489 
(u sin i -10 + 1 km s“'), there might be deviations from the con¬ 
stant frequency spacings. Fig. |9] indicates possible characteris¬ 
tic spacings in the range of 0.2-0.4, 1.0-1.4 and 2.0-2.7d“' in 
the low-frequency domain, and in the range of 0.4-0.7 and 2.6- 
2.7 d ' above 27 d ' in the high-frequency domain. 

We listed the possible rotational periods of KIC 9533489, 
calculated by Eq.|5]for some specific spacings in Table|5] These 
values correspond to local maxima of the ET of the frequencies 
below 10d“^ (green line in Eig.|9]) or above 27d '. We also con¬ 
sidered the possibility that in the case of a high inclination of the 
pulsation axis, or if any mode selection mechanism according to 
m is in operation, we may not detect the m - 0 components of 
the { - \ triplets, and the spacings belong to 6m — 2 splittings. 
Eor completeness, we also listed the possible rotational periods 
for i - 2 modes. We assume that the pulsation and rotation axes 
coincide. 

To check if the two dominant modes (/i and / 2 ) are rotation- 
ally split frequencies, we included their 0.44 d ' separation in 
Table|5] too. Finally, we included the case of /34 (1.5622d '). It 
takes part in two combinations :/39 = / 12+/34 and /45 = /i 2 -/ 34 - 
This practically means that there is an equidistant triplet struc¬ 
ture around /12 (6.3656 d '), with the separation of the value of 
/ 34 . There are no other similar structures found, but 1.56 d“^ is 
another frequency separation we have to take into account as the 
possible rotational spacing. 

According to Eq.|7l we expect that 6vp ^ 26vt=i^g. Consider¬ 
ing the 0.37 and 0.66 d“^ characteristic spacings of the presumed 
g- and p-modes, these might correspond to rotational splitting 
values of ^ = 1 modes. Another possibility is that the 1.30 and 
2.67 d“* spacings are caused by the stellar rotation. 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 


22 p- 
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Fig. 9. Searching for characteristic frequency spacings. From top to bot¬ 
tom: histograms of frequency differences, results of K-S tests, Fourier 
transforms of different frequency sets. Red boxes/lines: frequencies be¬ 
low 5 d“*, green boxes/lines: frequencies below lOd^*, blue boxes/lines: 
frequencies below 15 d“*, gray boxes/lines: frequencies above 27 d“*. 


None of these spacings selected in the gr-mode region shows 
clearly the 0.6 ratio expected for ( - I and 2 modes (see Eq.| 6 ll. 

Utilizing the rotation periods calculated from the different 
spacings, the estimated radius of the star, and the observed v sin i 
(see TableO, we indicated for guidance the corresponding equa¬ 
torial velocity (v) and inclination values in Table |5] This re¬ 
veals that since the lower limit for ti is a; 63kms“^ (vsini = 
70 - 7kms“', i - 90°, cf. TableO, the 0.37d“' separations of 
the presumed g-modes cannot be interpreted as rotational spac¬ 
ings. This also holds for the 0.66 d ' spacing of the presumed 
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Fig. 11. Searching for characteristic period spacings. Notations of the 
different period sets are the same as in Fig.|9] 


Fig. 10. Echelle diagrams of the 12 independent frequencies, supposed 
p-modes, above the frequency gap. 


p-modes, but not for the 0.44 d * frequency separation (cf. Ta- 
ble|5]l. Therefore 0.44d * is a plausible rotational spacing. 

The (critical) Keplerian velocity of tiKep = 437kms“* is 
reached at an inclination of i - 9.4°. It seems improbable that 
the 2.37 d“’ spacing represents rotationally split { - 1,6m = 1 
gr-modes, as the derived equatorial velocity of 371 km s ' is more 
than 80 per cent of tiKep- 

The 2.67 d * spacing found from the high-frequency re¬ 
gion requires further investigation. We have discussed already 
that it is a potential rotational spacing value. However, it can 
also be related to the large separation of frequencies. There 
are known cases in the literature, when such peaks in the 
Fourier transforms of the observed frequ encies turned out to 
be at one half of the large se paration (e.g. [Handler et al.lll99^ 
[Garcia Hernandez et al.ll2009l) . It is expected when m odes with 
different f values are observed. ISuarez et al.l (1201 4ll present a 
method to estimate the large separation (Av) for bSct stars 
(Eq. O. Applying their formula, in the case of KIC 9533489 
(p/po 0.4), the large separation is expected to be 68.6pHz 
(5.9d *). Half of this value is 2.95 d *, which is not far from 
the 2.67 d * spacing observed. Considering the uncertainties in 
the radius and mass determination, and that equidistant spacings 
are expected for slow rotators and frequencies in the asymptotic 
regime, the 2.67 d ' value is a reasonable estimate for half of the 
large separation. 

We demonstrate in Fig.[T0]that this 2.67 d * frequency spac¬ 
ing is not a mathematical artefact, but indeed exists in the fre¬ 
quency set. The plot shows echelle diagrams; the frequencies 
above the gap as a function of frequency modulo both the 
2.67 d“* spacing and its double. Forming vertical structures, 5-7 
frequencies may be associated with these spacings. 


In conclusion of this short discussion of the frequency spac¬ 
ings, we see two options here. One of them is that the pulsator 
has a high inclination (/ > 70°) and is a moderate rotator with an 
equatorial velocity of around 70kms *. In this case, frequency 
differences with » 0.4 d ’ in the g'-mode region can be related 
to rotational splitting. This also could explain the presence of 
the 0.84 d ' peak and its harmonics in the light curve’s Fourier 
spectrum, as being caused by rotational modulation. 

Another option is that KIC 9533489 is a fast rotator, with 
V » 200 kms * and i =3 20°. This could explain the triplet with 
1.56d ' frequency separation or the l.l-1.3d * characteristic 
spacings shown in Fig.j^as rotational spacings. In this case, the 
rotational spacing for the p-modes and half of their theoretical 
large separation would lie close to each other. 

5.2.3. Periocd spacings 

Using Eqs.[T]and|2] we performed similar tests on the periods 
as on the frequencies in Sect. 15.2.21 Only the histograms and 
the Fourier transforms suggest the presence of some characteris¬ 
tic spacings around 0.012-0.017, 0.002, 0.024-0.026 and 0.032- 
0.035 d (cf. Fig.fTTTi. We cannot use and interpret these indefinite 
results. Modelling of KIC 9533489 could help to decide if any of 
these spacing values indeed characterise { modes of the star. 

Note that period spacings are strongly influenced not only by 
chemical gradients, but also by the stellar rotation, especially in 
the case of modera te and fast rot a tors, a s demonstrated by the 
theoretical study of iBouabid et al.l (l2013h . 


5.2.4. Other regularities 

Coupling between g- and p-modes was re ported in many 
hybri ds. In the case of C oRoT 105733033 (IChape llier et al.] 
120121) . CoRoT 100866999 (IChapellier & MathiasI 120131). and 
KIC 11145123 (iKurtz et al.ll^l4 ). frequencies can be detected 
in the p-mode regime according to the relation fi p - F ± fi^g. 
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Table 5. Possible rotational periods (P), equatorial velocities (v) and inclinations (i) of KIC 9533489 applying different frequency spacing values. 


Low-frequency region (below 10 d 

Dipole modes (f = 1) Quadrupole modes (f = 2) 


Spacing 

P[d] 


V [km s ' ] 

i [deg] 

P[d] 


V [km s ] 

i [deg] 

[d-i] 

[pm 

6m—1 6m—2 

6m—1 

6m—2 

6m—1 

6m —2 

6m—l 6m—2 

6m=l 

6m—2 

6m—1 

6m—2 

0.37 

4.3 

1.4 

2.7 

57.9 

28.9 

- 

- 

2.3 

4.5 

34.7 

17.4 

- 

- 

1.30 

15.0 

0.4 

0.8 

203.4 

101.7 

20 1 

43.5^1 

0.6 

1.3 

122.0 

61.0 

35.0®;j 

- 

2.37 

21A 

0.2 

0.4 

370.7 

185.4 

10 9^^ *^ 

99 9 24.5 

0.4 

0.7 

222.4 

111.2 


39.0 

0.44 

5.1 

1.1 

2.3 

68.8 

34.4 

“66.2 

- 

1.9 

3.8 

41.3 

20.6 

- 

- 

1.56 18.1 0.3 0.6 244.0 

High-frequency region (above 27 d"') 

122.0 

16 7 

15.0 

35 0^^-^ 

0.5 

1.1 

146.4 

73.2 

28.6,3^.5^ 

23 . 059.4 

0.66 

7.7 

1.5 

3.0 

51.8 

25.9 

- 

- 







2.67 

30.9 

0.4 

0.7 

209.1 

104.6 

19 6^^-^ 

42 








Notes. Rotational periods were calculated according to Eq.|5] equatorial velocities are derived from R, and P, and inclinations from these v values 
and assuming v sin i = 10 ±7 km s“*. 


where f^g is a g-mode frequency and F is the dominant p-mode: 
the radial fundamental or a low-overtone frequency. The de¬ 
tection of such coupling supports the idea that the frequencies 
found in the low- and high-frequency region originate from the 
same pulsator. We note that KIC 11145123 is a remarkable ‘text¬ 
book example’ showing how the regular spacings and their de¬ 
viations can be used to study the stellar interior. A nother exam¬ 
ple of mode coupling is the case of KIC 8054146 (iBreger et alJ 
l 2012 h . where peaks were detected in the high-frequency region 
at /o + nfi, where /o is a constant and f\ is the dominant fre¬ 
quency in the low-frequency regime. 

We searched for such regularities in the frequency spectrum 
of KIC 9533489, but only one frequency in the p-mode regime 
can be explained as the combination of another high-frequency 
peak and a g-mode frequency, as /50 = /34 - 1 - /v (cf. Table nil, 
where /50 = 34.14, /34 = 1.56 and fj - 32.58 d"^ It is worth 
noting that /34 corresponds to the frequency separation of a 
triplet, as discussed in Sect. 15.2.21 

6. Transit events 

The light curve of KIC 9533489 shows transits, which are shal¬ 
low (0.5%), short (0.07 d) and recur every 197 days (it is also 
planetary candidate KOI-3783.01). 

First, we checked the field of KIC 9533489 searching for 
possible external sources, which could be responsible for the ob¬ 
served transit events. There are indeed two known stars offset by 
~ 4 and 8 arcsec that are blended within the Kepler photomet¬ 
ric aperture. However, these are both too faint to cause 0.5 per 
cent deep eclipses even if fully eclipsed. Additionally, the Kepler 
data centroid analysis for this 40cr transit event (averaged over 
all events) shows that the transit host is coincident with the target 
star within the 3(T limit of 0.25". Considering the adaptive optics 
observation (see Sect. Si, we know that there is an additional star 
ss 1" from the primary object. However, it appears to be too faint 
to be the transit host star, and it is also still too far away from the 
primary. 

We mode lled the Kepler L C and SC data using the jkte- 
Bop code (see lSouthworthll2008l and references therein), and ac¬ 
counted for the low sampling rate of the LC data by numerically 
integrating the model to match. Whilst the SC data cover only 
one transit, their better sampling rate makes them more valuable 
than the LC data, which cover six transits. We pre-whitened the 


data to remove the effects of the pulsations, then extracted the 
regions around each transit, and normalised them to unit flux by 
fitting a straight line to the out-of-transit data. 

The transits are indicative of a planetary-mass object orbiting 
the FO star, but we encountered problems with this hypothesis. 
The duration is far too short for this scenario and implies a much 
smaller radius for the FO star than found from our spectroscopic 
analysis. The low ratio of the radii implied by the transit depth 
also yields transit light curves with ingress and egress durations 
too short to match the data. These problems forced us to consider 
the possibility that the putative planet is orbiting a fainter star, 
which is spatially coincident with the FO star. This allows the 
use of a larger ratio of the radii, yielding a better match to the 
durations of the partial phases of the transit, and also implies a 
smaller stellar radius and thus a shorter transit (Fig. [13. 

Fitted models of the light curve including a larger ‘third 
light’ (L 3 ) component are ind eed significantly bet ter than those 
with L 2 - 0 (see Fig. [T3. ISouth^wo^ (12^012h showed that 
L 3 is very poorly constrained when modelling transits in iso¬ 
lation, but only considered values between 0.0 and 0.75 (ex¬ 
pressed as a fraction of the total light of the system). We find that 
L 3 = 0.88+0.03 yields a determinate and improved fit to the light 
curve. The fainter star should therefore be roughly eight times 
fainter than the FO star, which implies a mass of about 0.95 Mq 
if the two stars are at the same dista nce. We used the D SEP the¬ 
oretical stellar evolutionary models (iDotter et alJl2008l) to guide 
these inferences of the physical properties of the system. 

Whilst a 0.95 Mq star has a much smaller radius than the FO 
star, we still require a large orbital eccentricity for the transit 
model to match the observations: e sincu = 0.956 where e is or¬ 
bital eccentricity and o) is the argument of periastron. The large 
L 3 value we find means that a larger ratio of the radii (k = 0.185 
versus 0.064) is needed to match the transit depth. Such an ob¬ 
ject orbiting a 0.95 Mq main sequence star will have a radius of 
approximately Q.16Rq (L6/?jup). This is consistent with the ra - 
dius of an inflated transiting planet (see lSouthworth et al.ll2012ll . 
but also with that of a low-mass star. We therefore cannot claim 
that the transiting object is a planet. 

If we allow the transit host star to be closer to Earth than 
the FO star, making them merely an asterism rather than grav¬ 
itationally bound, then its contribution to the system light can 
be matched by a less massive and thus smaller star. This would 
mean that the light curve could be matched using a smaller 
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Fig. 12. Kepler LC (above) and SC (below) light curves compared to 
the best hts from jktebop neglecting third light. The residuals are plot¬ 
ted at the base of the figure, offset from unity. The blue and the purple 
lines through the LC data show the best-fitting model with and with¬ 
out numerical integration, and the ones through the residuals show the 
difference between this model with and without numerical integration. 
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Fig. 13. As Fig. [12] but allowing for third light when calculating the 
best ht. 

e sin w, and also imply a lower radius for the transiting object. 
We do not currently have sufficient observational constraints to 
constrain the relative distances of the two stars. Table |6]gives a 
set of plausible physical properties of the system, based on the 
transit light curves, the DSEP theoretical stellar models, the as¬ 
sumption that the FO star and the transit host star have a common 
distance, and that the transiting body is stellar rather than an in¬ 
flated exoplanet. 

7. Conclusions 

KIC 9533489 is one of the many A-F spectral type stars ob- 
se rved by Kepler. It was sele cted as a hybrid candidate pulsator 
bv lUvtterhoeven et al.l (1201 ih . What makes it a promising target 
for detailed analysis is its hybrid nature which is most probably 
genuine. Besides that, it also shows periodic transit-like events 
in its light curve. 

One of the relevant questions is whether KIC 9533489 is a 
genuine hybrid showing simultaneously self-excited p- and g- 
modes or not. The asteroseismic relevance of such genuine hy¬ 
brids is very high, as their structure can be probed down to 
the deep interior layers. Based on the spectroscopically derived 


Table 6. Plausible model of the KIC 9533489 system. 


Orbital parameters 

To (BJDtdb) 2455527.46780 + 0.00025 
P(d) 197.14468 + 0.00015 

esinw 0.956 

ecosm 0.0 fixed 


Photometric parameters 

ta = Rfi^la 

0.0060 fixed 

i (°) 89.97 

k - Rb/Ra 

0.185 

L 3 0.875 

tb = Rb/o 

0.0011 



Transit host star 

transiting body 

Mass (Mq) 

0.95 

0.15 

Radius (Rq) 

0.85 

0.16 

Teff (K) 

5500 

3200 


Notes. L 3 is a fraction of the total light in the Kepler passband which 
is not coming from the transit host star, a is the orbital semimajor axis 
and f?A,B are the true radii of the objects. 


physical parameters of KIC 9533489 showing its location within 
the overlapping region of the 6 Set and y Dor instability strips, 
on the dominant modes’ coupled amplitude variation, and on 
the strong multiperiodicity in the low frequency domain, we 
argue that KIC 9533489 is a genuine hybrid pulsator. Further¬ 
more, based on frequency spacing investigations, we come to the 
conclusion that KIC 9533489 is either a highly inclined mod¬ 
erate rotator (u » 70kms“', i > 70°) or a fast rotator with 
V ss 200 km s * and i « 20°. We need a comprehensive modelling 
of the star which, among other parameters, takes into account the 
effects of rotation. 

The transit analysis revealed that KIC 9533489 (or the 
KIC 9533489 system) is unusual. It seems, that we have to con¬ 
sider the possibility of a high-order system; the Kepler aperture 
contains (1) an FO star showing pulsations, (2) a fainter star at 1" 
(detected by the adaptive optics measurements), (3) an (unseen) 
G/K-type star in the AO PSF of the FO star (either gravitation¬ 
ally bound or unbound), and (4) a small/faint object transiting 
the G/K star. Since stellar binaries (or higher order systems) are 
common, it is certainly possible to have a hierarchical triple sys¬ 
tem (objects 1,3,4). Further radial velocity observations on long 
(years) time scale and higher signal-to-noise ratio spectra might 
help to define the spectral type of the unseen transit host star and 
thus help to characterise the transiting object, too. 
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